Late-Glacial and Early Holocene Environmental and Climatic Change at Lake Tambichozero, Southeastern Russian Karelia

High-resolution lithostratigraphy, mineral magnetic, carbon, pollen, and macrofossil analyses, and accelerator mass spectrometry ¹⁴C measurements were performed in the study of a sediment sequence from Lake Tambichozero, southeastern Russian Karelia, to reconstruct late-glacial and early Holocene aquatic and terrestrial environmental changes. The lake formed ca. 14,000 cal yr B.P. and the area around the lake was subsequently colonized by arctic plants, forming patches of pioneer communities surrounded by areas of exposed soil. A minor rise in lake productivity and the immigration of Betula pubescens occurred ca. 11,500 cal yr B.P. The rise in summer temperatures probably led to increased melting of remnant ice and enhanced erosion. The distinct increase in lake productivity and the development of open Betula-Populus forests, which are reconstructed based on plant macrofossil remains, indicate stable soils from 10,600 cal yr B.P. onward. Pinus and Picea probably became established ca. 9900 cal yr B.P.     

Oxygen fugacity and geochemical variations in the martian basalts: implications for martian basalt petrogenesis and the oxidation state of the upper mantle of Mars

The oxygen fugacity of the Dar al Gani 476 martian basalt is determined to be quartz-fayalite-magnetite (QFM) −2.3 ± 0.4 through analysis of olivine, low-Ca pyroxene, and Cr-spinel and is in good agreement with revised results from Fe-Ti oxides that yield QFM −2.5 ± 0.7. This estimate falls within the range of oxygen fugacity for the other martian basalts, QFM −3 to QFM −1. Oxygen fugacity in martian basalts correlates with ⁸⁷Sr/⁸⁶Sr, ¹⁴³Nd/¹⁴⁴Nd, and La/Yb ratios, indicating that the mantle source of the basalts is reduced and that assimilation of crust-like material controls the oxygen fugacity. This allows constraints to be placed on the oxidation state of the martian mantle and on the nature of assimilated crustal material. The assimilated material may be the product of early and extensive hydrothermal alteration of the martian crust, or it may be amphibole- or phlogopite-bearing basaltic rock within the crust. In either case, water may play a significant role in the oxidation of basaltic magmas on Mars, although it may be secondary to assimilation of ferric iron-rich material.     

Methodology for construction, calibration and validation of a national hydrological model for Denmark

An integrated groundwater/surface water hydrological model with a 1 km² grid has been constructed for Denmark covering 43,000 km². The model is composed of a relatively simple root zone component for estimating the net precipitation, a comprehensive three-dimensional groundwater component for estimating recharge to and hydraulic heads in different geological layers, and a river component for streamflow routing and calculating stream–aquifer interaction. The model was constructed on the basis of the MIKE SHE code and by utilising comprehensive national databases on geology, soil, topography, river systems, climate and hydrology. The present paper describes the modelling process for the 7330 km² island of Sjælland with emphasis on the problems experienced in combining the classical paradigms of groundwater modelling, such as inverse modelling of steady-state conditions, and catchment modelling, focussing on dynamic conditions and discharge simulation. Three model versions with different assumptions on input data and parameter values were required until the performance of the final, according to pre-defined accuracy criteria, model was evaluated as being satisfactory. The paper highlights the methodological issues related to establishment of performance criteria, parameterisation and assessment of parameter values from field data, calibration and validation test schemes. Most of the parameter values were assessed directly from field data, while about 10 ‘free’ parameters were subject to calibration using a combination of inverse steady-state groundwater modelling and manual trial-and-error dynamic groundwater/surface water modelling. Emphasising the importance of tests against independent data, the validation schemes included combinations of split-sample tests (another period) and proxy-basin tests (another area).     

Stochastic Structural Modeling

A consistent stochastic model for faults and horizons is described. The faults are represented as a parametric invertible deformation operator. The faults may truncate each other. The horizons are modeled as correlated Gaussian fields and are represented in a grid. Petrophysical variables may be modeled in a reservoir before faulting in order to describe the juxtaposition effect of the faulting. It is possible to condition the realization on petrophysics, horizons, and fault plane observations in wells in addition to seismic data. The transmissibility in the fault plane may also be included in the model. Four different methods to integrate the fault and horizon models in a common model is described. The method is illustrated on an example from a real petroleum field with 18 interpreted faults that are handled stochastically.     

Initiation of the last glaciation in northern europe

The bio- and chronostratigraphy of the Eemian interglacial (marine isotope substage 5e) and an Early Weichselian glaciation (5d-a) established from representative and detailed sequences can be correlated with the deep-sea oxygen isotope stratigraphy, ice-core data, sea-level fluctuations and coupled ice sheet-climate models. Biostratigraphic sequences from Fennoscandian key sections are correlated with reference sequences from Estonia and from sections located near or beyond the margins of the last glaciation. Organic sediments previously attributed to Early and Middle Weichselian interstadial periods in Finland are argued to be redeposited and mixed older (last interglacial) material. Pollen and diatom spectra of the undisturbed materials suggest that the Eemian climatic optimum was followed by a continuously cooling climate and a regressive marine level. If only undisturbed sequences are considered, the major climatic fluctuations of the Early Weichselian, apparent in Central and Western Europe, are not apparent in the sequences from the central part of the glaciated terrain. Instead, some sequences are truncated by sediments indicating approaching ice sheets soon after the interglacial. This may imply that the ice sheet grew over Finland during the first Early Weichselian stadial. The preservation of the interglacial beds and the lack of younger non-glacial sediments support the interpretation that the area remained ice-covered until the final deglaciation. During the Early Weichselian, the Norwegian coast was probably occasionally ice free, similar to the coastal zone of Greenland today. The authors' interpretation of the Fennoscandian organic deposits of the last glaciation may also explain similar observations from the central parts of the Laurentide ice sheet.     

Using rainfall-runoff modeling to interpret lake level data

Using water balance computations, the behavior of different kinds of lakes is discussed. Simple analytical expressions relating water level to hydrological conditions and lake bathymetry are given. The importance of knowing the river basin area when analyzing lake levels is stressed. A conceptual rainfall-runoff model including lake routing is used to simulate runoff and lake levels and to compute quasi-steady state conditions and long-term transient situations. It is suggested that models can be used to construct curves relating lake levels to precipitation and lake evaporation. By comparing with paleo-lake levels, the annual precipitation related to these levels can be found, provided information is available about the seasonal distribution of the precipitation.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Impact of warmer climate on melt and evaporation for the rainfed, snowfed and glacierfed basins in the Himalayan region

The impact of warmer climate on melt and evaporation was studied for rainfed, snowfed and glacierfed basins located in the western Himalayan region. Hydrological processes were simulated under current climatic conditions using a conceptual hydrological model, which accounts for the rainfall–runoff, evaporation losses, snow and glacier melt. After simulations of daily observed streamflow (R²=0.90) for 6 years, the model was used to study the impact of warmer climate on melt and evaporation. Based on the future projected climatic scenarios in the study region, three temperature scenarios (T+1, T+2 and T+3 °C) were adopted for quantifying the effect of warmer climate. The comparison of the effect of warmer climate on different types of basins indicated that the increase in evaporation was the maximum for snowfed basins. For a T+2 °C scenario, the annual evaporation for the rainfed basins increased by about 12%, whereas for the snowfed basins it increased by about 24%. The high increase of the evaporation losses would reduce the runoff. It was found that under a warmer climate, melt was reduced from snowfed basins, but increased from glacierfed basins. For a T+2 °C scenario, annual melt was reduced by about 18% for the studied snowfed basin, while it increased by about 33% for the glacierfed basin. Thus, impact of warmer climate on the melt from the snowfed and glacierfed basins was opposite to each other. The study suggests that out of three types of basins, snowfed basins are more sensitive in terms of reduction in water availability due to a compound effect of increase in evaporation and decrease in melt. For a complex type of basin, the decrease in melt from seasonal snow may be counterbalanced by increase in melt from glaciers. However, on long-term basis, when the areal extent of glaciers will decrease due to higher melt rate, the water availability from the complex basins will be reduced.